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Abstract 

Starting from the the Boltzmann equation, we study the center of mass oscillation of a harmon- 
ically trapped normal Fermi gas in the presence of a one-dimensional periodic potential. We show 
that for values of the the Fermi energy above the first Bloch band the center of mass motion is 
strongly damped in the collisional regime due to umklapp processes. This should be contrasted with 
the behaviour of a superfluid where one instead expects the occurrence of persistent Josephson-like 
oscillations. 
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It is well known that the center of mass of an atomic cloud confined in a harmonic potential 
oscillates with the trap frequency irrespective of temperature, inter-particle interactions and 
quantum statistics (Khon's theorem). This general result is no longer true if the system is 
also confined by a periodic potential generating barriers separating different wells. While 
a Bose-Einstein condensate still exhibits a collective behaviour since atoms can tunnel in a 
coherent way through the barriers Q], a non interacting Fermi gas with Fermi energy lying 
above the first Bloch band is not able to oscillate around the minimum of the harmonic well 
Pi- The sudden shift of the harmonic trap indeed causes an asymmetry in the occupation 
of the left and right orbits which, due to their open nature, persists in time. This causes the 
system to remain trapped at one side of the harmonic field. The inclusion of interactions 
changes this scenario favouring the relaxation towards equilibrium P|. 

An important question is to understand what is the behaviour of a normal Fermi gas 
in the coUisional regime where the system is in conditions of local equilibrium and the 
dynamic behaviour is described by the equations of hydrodynamics. In particular the main 
question addressed in the present letter is whether in the collisional regime the system is 
able to exhibit center of mass oscillations. The problem is relevant in view of the search of 
signatures of superfiuidity in ultra cold Fermi gases. 

We consider a two-component atomic gas trapped by an external potential given by the 
sum of a harmonic trap of magnetic origin and of a stationary optical potential modulated 
along the 2;-axis. The resulting potential is given by 

Vext = {ujl_rl_ + ulz^) + sEr siv? qbz, (1) 

where uj±, uJz are the frequencies of the harmonic trap, Er = /j,^g^/2m is the recoil energy, 
qb being the Bragg momentum and s is a dimensionless parameter providing the intensity 
of the laser beam. The optical potential has periodicity d = Ti/qs along the z-axis. In the 
semiclassical approximation, the energy of the resulting Bloch states is given by e(r, p) = 
e{pz) +p\/1m + V/io(r), where Vho = m{ijj']r\ + u;^2;^)/2 and e{pz) is the dispersion relation 
of the lowest band obtained by solving the ID Schrodinger equation with the Hamiltonian 
H = pl/2m + sERsirL^^qsz). Since e{pz) is a periodic function with period 2hqB, the quasi- 
momentum pz is restricted to the first Brillouin zone defined by —hqB < Pz ^ ^qs- We 
restrict our discussion to low temperatures T <^Tp and neglect higher bands. We consider 
situations where the system is in the collisional (hydrodynamic) regime and we calculate the 
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relaxation rate of the center of mass oscillation due to umklapp collisions. In order to achieve 
the hydrodynamic condition at such low temperatures, where collisions are suppressed by 
the Pauli principle, one should increase the value of the scattering length working close to a 
Feshbach resonance and/or work with very shallow traps along the z-th direction. 

We begin our analysis by defining the center of mass coordinates and momenta as Z = 
J zfdpdr/h^, Pz = J Pzfdpdr/h^ where / = f^ + fi and /|, fi are the distributions functions 
of the two spin species which are assumed to be equal (/^ = = //2). 

By suitable integrations of the Boltzmann equation, one finds the following exact equa- 
tions for the center of mass oscillation: 

—Pz + mUzZ = J PzC-^ , (3) 

where C is the coUisional integral describing s— wave scattering between Bloch states. Eq. Q 
directly follows from the equation of continuity, jz = de/dpz being the current density along 
z. In the absence of the periodic potential, the integral in Eq.Q is zero as momentum is 
rigorously conserved during collisions. Since in this case the dispersion law reduces to the 
free value p1/2m, one has J de/dpzfdpdr/h^ = Pz/m and one recovers the general result 
uj = Uz for the frequency of the oscillation. 

The effect of the optical lattice on the above equations is twofold. First it changes the 
dispersion relation from the free value p1/2m to e(p^). Second, the traslational symmetry 
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is broken and momentum is no longer a good quantum number 
consequences for the collisional integral as momentum conservation is replaced by the weaker 
constraint piz + P2z — Psz — Paz = 2hnqB, where n is an integer and piz,P2z and P3z,P4,z are, 
respectively, the initial and final quasi-momenta of the two colliding particles. For two-body 
interactions, one finds that only umklapp processes with n = ±1 are allowed. 

In an umklapp collision, the system exchanges momentum ±2hqB with the optical lattice. 
This means that Pz varies in time not only because of the oscillator force F = —muj'^z but 
also as a result of umklapp processes, which act to damp the oscillation. We therefore expect 
that the right hand side of Eq.® is non-zero for out-of-equilibrium distributions. 

Let us write the collisional integral as C(p, r) = C+(p, r) — C~(p, r), where C"'"(p, r) and 
C^(p, r) describes collisions in which one of the particles has, respectively, final and initial 
momentum p. Introducing the compact notation dp = dp±Q{\hqB\ —Pz)dpz, where Q{x) is 



the step function, the term can be conveniently expressed via the Fermi golden rule in 
the general form C+(pi, r) = J dp2dp3dp4D , with 
n( . 2 27C 

\Uj^^\'^6{piz+P2z-P3z-P4z-'2nqB) 

n=0,±l 

5(p±i + P±2 - P±3 - P±4)5(e(l) + e(2) - e(3) - e(4)) 
/T(3)/i(4)(l-/T(l))(l-/i(2)), (4) 

where Uj^^ is the n— dependent matrix element of the two-body interaction and e(j) = 
e(r, pj). The factor 2 in Eq.Q comes from the trace over spin indices and fa are normalized 
to J ffjdpdr/h^ = N/2. The C~ term is simply obtained from by interchanging initial 
and final states. At thermal equilibrium one has fa = fo = [e^^^~^^ + 1]^-*^ with /? = l/ksT 
and E = e(r,p). In this case = C~ and hence C = 0. The total number of collisions 
per unit time is given by 

^ ~ j J/z ~ ^noT + ^uk} (5) 

where we have written explicitly the contributions Tnor and Tuk coming from normal {n = 0) 
and umklapp {n = ±1) collisions [see Eq.(jH)]. 

While in the general case Eqs. and (jHl) are not sufficient to calculate the frequency of 
the mode, in the hydrodynamic regime one can make the ansatz 

/a(r, p, t) = /o(e(r, p) + u{t)p,), (6) 

where u{t) is a time dependent function. This corresponds to a rigid displacement of the 
density current jz'- jz{Pz) — ^ jz{Pz) + u and permits to write Eqs and (jHl) in a closed 
form. 

The ansatz (jH)) is valid in the limit of strong collisions uor ^ 1, where uj is the frequency 
of the collective oscillation and r is a typical coUisional time. In the classical regime one has 
T /N . At low temperatures T -C Tp, only the particles near the Fermi surface can be 
scattered, because of Fermi statistics, and one finds that the hydrodynamic condition takes 
the form 

- ~ > uj. 7 

Notice that at low T the rate F behaves like (T/Tp) [5] so that 1/r ~ [T/Tp) coincides, 
apart from a numerical factor of the order of unity, with the quasiparticle lifetime calculated 
at the Fermi surface 



Since both the integrated current J j^dz and the colhsional integral are zero at equihb- 
rium, we expand the ansatz ® to first order in u and plug the result into Eqs.Q and Q. 
Analogously, for the center-of-mass momentum one finds Pz = J Pzfdpdr/h^ = urfiN, where 

2 f 2 9fo dpdr 



plays the role of an effective mass for the center of mass oscillation |8[. This permits us to 
cast Eqs.(j21) © in the following closed form, holding in the limit of small amplitude 
oscillations: 

(9 P 

ot m 

^R + mu;lZ=-^, (10) 

where 

1 1 / PizjPiz + P2z - P3z - P4z)Ddpidr 

Tuk ksT J pldfo/de{pz)dpidr 

defines the relevant relaxation time of the oscillation due to umklapp collisions. To derive 
Eq-dlH) we have used the identity J2'j=iPj(^P' /^^iPjz) = -Pipiz +P2z-P3z-PAz)P>, holding 
at equilibrium also for Fermi statistics. Notice that the contribution from the normal {n = 0) 
collisions, which conserve momentum, identically vanishes. Equations and f lTUI) have the 
form of a damped harmonic oscillator. Looking for solutions of the form e"""^* the dispersion 
law is given by 

u = —±}-ul ^ (12) 

showing that the oscillations become overdamped if oJ^Tuk < 

By recalling that the integrand D in Eq.(0]) is symmetric under interchange 1^2 and 
antysimmetric under interchanges 1 ^ 3 or 1 ^ 4, the factor pi^ in the numerator can 
be substituted by the combination {piz + p2z — Psz — Piz)/'^ which is equal to ±hqB/'2 for 
(n = ±1) umklapp processes. Taking into account Eq.®, Eq. fjllj) can then be written as 

1 _m En Vuk 

— =A—— — (13) 
Tuk mknT N 

where F^^ is defined by Eqs.dH) and (jS}- 

Our next goal is to show that, for sufficiently tight optical lattices and gas densities 
corresponding to Tp/26 > 1, the condition ((Tj) automatically implies the overdamping of 
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the center of mass oscillation. In the following we consider relatively large laser intensities 
and work in tight-binding approximation. Under this assumption, the energy dispersion is 
given by e{pz) = 6(l — cos[^]) where the bandwidth 26 is proportional to the tunneling 
rate between consecutive wells. The Fermi energy is related to the total number of particles 
by = 2 / Q{Tf - e(p, r))dpdr/h^, yielding 

where p = Pzd/h and h{p) = Tp/6 — 1 + cos p. Eq. p4|) permits to calculate Tp as a function 
of the free Fermi energy Tp = {3N)^^^h{uj±'^uJzy^^ , evaluated in the absence of the optical 
lattice, the recoil energy En and the bandwidth 26. Neglecting higher bands, the Fermi 
energy Tp is always smaller than the free value Tp reflecting the fact that e{pz) < pi/ 2m in 
the first Brillouin zone. 

In order to calculate the relevant interaction matrix elements f//j, we first neglect har- 
monic trapping. The scattering states are e'^'^^''^^ipk{z), where ipk{z) is the ID Bloch wave- 
function with quasi- momentum k. We make the ansatz ipk{z) ~ 

/ labels the wells, / is localized at the origin and normalized to J^j^ \f{z)\'^ = l. In the 
following we will consider a delta function potential U{r) = gS{r), the coupling constant g 
being related to the scattering length a hj g = Anh^a/m. In the tight-binding limit, the 
matrix elements are given by 

^/f ^ = = 9d / f{z)dz = ga. (15) 

J-d/2 

Conversely, in the absence of the optical lattice the eigenstates of H are simply plane waves 
and the matrix elements are given by f^jf^"* = and f/j°^ = g. The factor a of Eq. ()15|l 
is larger than one and increases by increasing the laser intensity s, the wavefunction f{z) 
becoming more and more peaked at the origin. For very large values of s, one can use the 
asymptotic result a = ^J'n I2s^l'^ . This formula actually overestimates the correct value of 
a for smaller s. As an example, for s = 5 it gives a = 1.9 while a more accurate calculation 
gives a = 1.6 7]. From the above discussion, we conclude that the periodic potential 
enhances interaction effects in two different ways. First it allows for additional (umklapp) 
collisions to take place. Second, the optical confinement compresses the gas and this results 
in an increase of the coupling constant from g to ag. 

In the presence of harmonic trapping, we can still use result p5|) provided local density 
approximation is applicable. This requires that the trapping frequencies should be small 
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compared to the Fermi energy and the bandwith. The coUisional rates Tnor and Tuk ap- 
pearing in Eqs. (0} and (0) have been integrated using standard Montecarlo techniques. The 
ratio Tuk/Tnor is plotted in Fig^as a function of the temperature for different values of the 
parameter Tp/26. For Tp <^ 26, umklapp collisions are negligible at low temperatures. In 
fact, in this case, the typical initial quasi-momenta piz,P2z are small compared to hqB and, 
due to the energy conservation, this is also true for the final quasimomenta, meaning that 
processes with n are unlike. When the Fermi energy is larger than the bandwith, umk- 
lapp processes become instead competitive with the normal ones even at low temperatures 
and we see that the ratio Tuk nor saturates to a constant value. In a uniform system this 
constant can be analytically evaluated and is equal to 1/2. In the trapped case one finds a 
smaller value because the effective Fermi energy Tp{r) = Tp — Vho{r) is r— dependent and 
therefore umklapp collisions are quenched at the edge of the cloud where Tp{r) < 26. 

We are now ready to compare the overdamping condition with the hydrodynamic con- 
dition ((Zj). To this purpose, we have evaluated the effective mass m for the center of mass 
oscillation introduced in Eq.(jHl). At T = and for a fixed laser intensity s, the ratio m/m is 
a function of the parameter Tp/En. This function is plotted in Fig.Q for different values 
of s. The figure shows that, for Tp/Eji > 1, the ratio m/m does not depends on the laser 
intensity and is of the order of unity. 

In the following we will consider typical configurations with Tp ~ Ef{. By comparing 
Eq.((Zj) with Eq. (fT3|) . we conclude that, due to umklapp processes, the center of mass oscilla- 
tion of a trapped gas confined by a tight optical potential with Tp > 26 will be overdamped 
in the collisional regime cj^r -C 1 since in this case Tuk ~ Tnor, m m and therefore 
Tuk ~ r < u;^^ 

The relaxation rate of the dipole oscillation can be conveniently written in the form: 



where the temperature dependence of the function F is plotted in FigOlfor different values 
of the paramter Tp/26. The possibility for the system to be in the overdamping regime 
depends on the actual values of the parameters in Eq. (fT6|l . As a concrete example, let us 
consider a two-component gas of = 10^ potassium {'^^K) atoms with trap frequencies 
uj± = 27T ■ 275Hz and uj^ = 27t ■ 24Hz, corresponding to Tp = 390nK. For the optical 
lattice we assume s = 5 and periodicity d = 400nm, corresponding to Er = 9.26 = STinK. 




(16) 



7 



For the Fermi energy one then finds Tp = 0.85Tp and from the value of Tf/2S one finds 
m/m = 1.8. We see from FigjSlthat the overdamping condition oOz^uk <^ 1 is satisfied even 
at low temperatures T/Tp ~ 0.05 — 0.1, provided the scattering length is moderately large, 
say \a\/d > 0.1 — 0.2. This can be accomplished experimentally working close to a Feshbach 
resonance. 

For lower temperatures or smaller scattering lengths, the system enters the coUisionless 
regime. Since in this regime the center-of-mass oscillation is self-trapped under the same 
Tp > 26 condition [^, we conclude that the system can never exhibit undamped center-of- 
mass oscillation in the normal phase. In the superfiuid phase, on the contrary, one expects 
the occurrence of persistent (undamped) Josephson-like oscillations . In particular, in the 
weak-coupling (BCS) limit, where the distribution function / does not deviate significantly 
from the ideal gas value, the frequency of the dipole oscillation is expected to be ujz\fmf^ 
with m given by Eq.® and Fig|21 

Let us finally discuss the conditions of applicability of our approach. These concern 
the stability of the Bloch states with respect to the interaction. First, the scattering length 
must be small compared to the interwell distance, i.e. a <^ c/, otherwise we cannot model the 
interaction with a 5-function potential. Second, the broadening of the Bloch wavefunction 
due to collisions should be smaller than the bandwith: h/r <^ 5. This condition is needed 
in order to apply the semiclassical Boltzmann picture and can be reasonably well satisfied 
for the temperatures of interest. 
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FIG. 1: Ratio of umklapp over normal collisional rates [see Eq.@] as a function t = T/Tp for 
the values Tp/26 = 5 (solid line), 2 (dashed line), 1 (dotted line). 
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FIG. 2: Effective mass of the dipole oscillation [see Eq.Q] plotted at T = as a function of the 
parameter Tp/Epi for laser intensities s = 5 (solid line) and s = 8 (dashed line). The asymptotic 
curve bEfi/STp is also shown (dotted line). 
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FIG. 3: Dipole relaxation function F [see Eq. (|TH|) ] as a function of the reduced temperature 
t = T/Tp for the values Tp/25 = 5 (solid line), 2 (dashed line), 1 (dotted line). 
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